Non- Fermi liquid properties of 2d symplectic fermions: the role of 
a dynamically generated (pseudo)-gap 

Eliot Kapit and Andre LeClair 
Newman Laboratory, Cornell University, Ithaca, NY 

Abstract 

The interacting symplectic fermion model in two spatial dimensions is further analyzed. As an 
effective low energy theory, the model is unitary. We show that a relativistic mass m is dynamically 
generated and derive a gap equation for it. By incorporating a finite temperature we study some 
fundamental properties of the model, such as the specific heat and spin response, which clearly 
show non-Fermi liquid properties. We find that various physical properties are suppressed at tem- 
peratures T < T* where the cross-over scale is T* = m. As a simplified, toy model of high Tc 
superconductivity, we thus identify the pseudogap energy scale with the zero temperature relativis- 
tic mass m, and show that this reproduces some qualitative aspects of the observed phenomenology 
of the pseudogap. The effects of the pseudogap and finite temperature on the d-wave gap equation 
are analyzed. In this model, the pseudogap is a distinct phenomenon from superconductivity and 
in fact competes with it. Our analysis of Tc suggests that the quantum critical point of our model, 
where the pseudogap vanishes, occurs inside the superconducting dome near optimal doping. For 
an antiferromagnetic exchange energy oi J /ks ^ 1350i^, solutions of the d-wave gap equation give 
a maximum Tc of about 
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I. INTRODUCTION 



An important aspect of strongly-correlated electron physics is the expectation that some 
systems may exhibit novel non-Fermi liquid behavior. This can present some interesting 
challenges for theoretical models since non-Fermi liquid behavior is known to be rare based 
on simple renormalization group (RG) arguments[l|. The best understood example is the 
Luttinger liquid in one spatial dimension, where the non-Fermi liquid behavior is essentially 
attributed to the fact that quartic interactions of Dirac fermions are marginal in the RG 
sense. In higher dimensions quartic interactions of Dirac fermions are irrelevant. These 
considerations were one of the primary initial motivations for the construction and anal- 
ysis of a new model of interacting fermionsj^, [s], with non-Fermi liquid behavior. Like 
the Luttinger liquid, the model is a continuum field theory of only 4 fermionic fields xfi 
which carry charge and spin, with a unique quartic interaction due to Fermi statistics. The 
novelty of the model is the free kinetic term in the hamiltonian which in addition the com- 
monplace term that is second order in spatial derivatives Vx~ ■ Vx^, contains an additional 
contribution that is second order in time derivatives dtX~dtX^ i ^i-nd thus has an emergent 
Lorentz symmetry. The model can be consistently canonically quantized as a fermion, even 
though this structure of the kinetic term is usually associated with relativistic bosons. In 
two spatial dimensions the field x has dimension 1/2, which implies the quartic interaction 
has dimension 2 and is thus RG relevant. The model in fact has a low-energy fixed point, 
i.e. quantum critical point. 

Several interesting properties in addition to the non-Fermi liquid behavior emerged in 
the analysis of the modelj4|. It has a hidden 5*0 (5) symmetry which contains the spin 
5*0(3) and electric U{1) as commuting subgroups. The fundamental fields i transform in 
the 4-dimensional spinor representation of 5*0(5) and the bilinears decompose as 4 ® 4 = 
1 © 5 © 10. The 5 vector representation serves as order parameters for the spontaneous 
symmetry breaking of 50(3), i.e. magnetic order, and also contains Cooper pair fields of 
charge ±2 which are order parameters for symmetry breaking of U (1), i.e. superconductivity. 
By deriving separate gap equations and studying their solutions it was shown that the model 
can contain an anti- ferromagnetic phase (AF) and a d-wave superconducting phase (SC). In 
this model the basic mechanism that leads to a d-wave SC instability is clearly identified as 
arising from the 1-loop scattering of Cooper pairs. An attractive feature of the model is its 
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simplicity, reflected in the fact that it has very few free parameters: an overall energy scale Eq 
set by the high-energy cutoff Ac, a Fermi velocity vp, and a single coupling g. Furthermore, 
at low energies g is near the low-energy flxed point value of 1/8, so that perturbation theory 
can be reliable. 

The above features led us to propose this theory as a toy model of high temperature super- 
conductivity (HTSC). Here the 5*0(5) symmetry is quite different from the S0{5) symmetry 
of Zhang in this context [s]: whereas he postulates an 5*0(5) invariant Landau- Ginsburg the- 
ory for the bosonic spin and electric order parameters, our model is a microscopic model of 
fermions where the bosonic order parameters are composite bilinears in the fermions. The 
5*0(5) symmetry was not put in by hand in our model, rather it is accidental and hidden. In 
fact, the AF and d-wave SC phases are not related by the 50(5) in our model. In principle 
a Landau- Ginsburg effective theory can be derived from our model, but like the BCS the- 
ory, our underlying microscopic theory is more powerful for deriving gap equations, critical 
temperatures, etc. Arguments that motivate its application to HTSC were given in ^|, the 
best one being its relation to the non-linear 0(3) sigma model effective theory for the 2d 
Heisenberg model description of the AF phase, which is also a relativistic theory; however 
a rigorous derivation of it from say the Hubbard model on the lattice is lacking. If our 
model really turns out to correctly capture some essential features of HTSC, then it should 
be viewed as an effective low-energy theory for wavelengths that are long compared to the 
lattice spacing. One expects that such a description can be approximately described by 
a rotationally invariant quantum fleld theory describing the gas of quasi-particles that SC 
condenses out of. The shortcoming of such a theory is that lattice effects, such as features 
that are dependent on the detailed structure of the Fermi surface in the Brillouin zone, are 
necessarily absent in the basic model, although perhaps they can be incorporated with small 
perturbations of the hamiltonian. Furthermore, our model does not have an explicit Fermi 
surface, as for the 0(3) sigma model of the AF phase at half-flUing. Nevertheless, the exis- 
tence of a simple model of fermions with only repulsive quartic interactions that possesses 
the all of the main motifs of HTSC, which is to say it has clearly identifled mechanisms for 
the d-wave SC, the pseudogap, and anti-ferromagnetism, can be useful for deciding which 
properties are essential and which are superfluous for the phenomenon. Thus, although it 
has not been established that our model captures what is really happening in HTSC, for 
the remainder of this paper we will take the liberty to borrow the terminology of the HTSC 
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literature where appropriate. 

The main purpose of the present article is study in detail the origin and properties of 
the so-called pseudogap in our model. In the HTSC materials, the pseudogap refers to an 
energy scale Epg = ksT* where a cross-over behavior is observed in a variety of physical 
properties such as electronic specific heat, magnetic susceptibility and conductivity. The 
same energy scale can be observed as the onset of a depression in the density of states. In 
the underdoped region, T* is considerably larger than the superconducting Tc. For reviews 
see Is, 7]. The origin and physical interpretation of the pseudogap remains a fundamental 
question in the physics of HTSC, and many researchers feel that a proper understanding of 
it will be an important key toward unraveling the mysteries of HTSC. Two broad classes 
of theories can be summarized by the schematic phase diagrams shown in Figure [H In the 
scenario on the left, T* goes to zero inside the SC dome, perhaps terminating at a quantum 
critical point. In this class of theories the pseudogap is unrelated, and in fact competes 
with, superconductivity. In the second class of theories, the pseudogap indicates pre-formed 
Cooper pairs for example, i.e. is a friendly precursor to SC, and the pseudogap line merges 
with Tc on the overdoped side. These opposing scenarios are discussed in some detail in 
js]. Thermodynamic data such as specific heat favors the first scenario [o], |lo|, whereas 
spectroscopic data seems to support the second 
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FIG. 1: Two proposed theoretical scenarios for the HTSC phase diagram. 

In our model there is an obvious identification of the pseudogap: whereas AF or SC order 
are related to the order parameters in the 5 vector of 5*0(5), the pseudogap is naturally 
associated with the singlet bilinear in the tensor product of 4 (g) 4, as suggested by Tyejl2|. 
and also more tentatively in jj]. This corresponds to the operator x~X~^ = J2a=-\ iXaXt- 
Whereas this term is classically zero in the relation with the 0(3) sigma model, it is dy- 
namically generated in the presence of interactions. As a contribution to the hamiltonian, it 
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corresponds to a relativistic mass m, and the 1-particle states have energy Ei^ = y/W+lv^. 
Thus a non-zero mass m gives a clean gap in the density of states, i.e. there is nothing 
"pseudo" about it. As discussed above, the anisotropy of the observed pseudogap is not a 
feature of this mass-gap since it is a lattice effect ("Fermi arcs"). In the sequel we identify 
this mass with the pseudogap energy scale Epg = T* = m and will show that its properties 
closely parallel the phenomenology of the pseudogap in HTSC, at least on the underdoped 
side. 

The pseudogap has an entirely different origin than the SC gap in our model, essentially 
because it corresponds to a dynamically generated vacuum expectation value (x'X^) which 
preserves the 5*0 (5) symmetry whereas SC breaks the electric U{1) subgroup. It is an 
intrinsic property of the normal state density of states. Whether it competes with or aids 
superconductivity is straightforwardly addressed by incorporating it into the d-wave gap 
equation derived in [4]. As we will see, it clearly competes with superconductivity since if 
it is too large it destroys the solution to the SC gap equation. Fortunately, the pseudogap 
goes to zero at the critical point and this makes superconductivity possible; in fact SC is 
enhanced near the critical point. Thus our model is in the class of the left figure in Figure [H 
with T* terminating at a quantum critical point near optimal doping. (The proper treatment 
of temperature presented in this paper led to the identification of the critical point which 



differs from the original proposal in 



m 



Our analysis of the pseudogap required two technical improvements of the work presented 
^. In the latter work, temperature was treated crudely as a mass, the idea being 
that temperature and a mass can both be viewed as an infra-red cutoff; this is ultimately 
unreliable for the computation of thermodynamic properties. In the present work, since the 
mass is identified with the pseudogap which is a zero temperature property, temperature 
must be dealt with properly, and this is accomplished with the Matsubara formalism. A 
non-zero mass also resolves some infra-red divergences that were present in the previous 
treatment. 

The remainder of this article is organized as follows. In the next section we review the 
definition of the model and some of its basic properties. Section III addresses the unitarity 
of the model, beyond what is contained in our previous work. Section IV describes our RG 
prescriptions. In section V we analyze the dynamical generation of a mass, i.e. pseudogap, 
by deriving a gap equation for the vacuum expectation value {x^X^)- The pseudogap 
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depends on the variable x which up to a scale is the inverse coupling. In section VI we 
propose a relation between the variable x and hole doping which thus gives the doping 
dependence of the pseudogap. In section VII we compute the effect of the pseudogap on 
the electronic specific heat. A non-zero magnetic field is introduced in section VIII, and we 
compute the temperature dependent spin susceptibility and the magnetic field dependence 
of the specific heat. For all of these thermodynamic properties, we find crossover behavior 
at the temperature T = T* and its qualitative dependence on doping compares favorably 
with data. Finally in section IX we derive the finite temperature version of the d-wave 
gap equation and incorporate the effect of the pseudogap into it. Analysis of this equation 
clearly shows that the pseudogap competes with SC, and leads to a computation of the 
superconducting T^. 



II. REVIEW OF THE MODEL, IT'S SYMMETRIES AND ORDER PARAME- 
TERS 

As for any second-quantized description of electrons with spin i, the fundamental fields 
of the model are 4 fermionic fields x^i where the flavor index a =t, i corresponds to spin 
and ± is electric charge. Due to the fermionic statistics there is a unique quartic interaction, 
thus various models are primarily characterized by the free kinetic term. Our model in two 
spatial dimenions is 



F = / d^x I ^ {dtX-dtXl + v'pVx- ■ ^Xl + m'XaXl) + Svr^t? Xf Xf Xl^i ) 

The above hamiltonian would be a standard second-quantized field theory for fermions 
interacting via a delta-function potential if it weren't for the term that is second order in 
time derivatives, and this is the primary novelty of the model. This choice of kinetic term can 
be motivated from the phenomenology of HTSC, since it leads to the correct temperature 
dependence of the specific heat C oc at low temperatures (see sectin VII) in the abscence 
of superconductivity, which is characteristic of a relativistic theory, and the mass m can 
correspond to the pseudogap. 

As a model of HTSC, the above kinetic term can also be motivated as follows^. Sup- 
pose one is near the Mott-Hubbard insulating phase. The anti-ferromagnetic phase of the 
Heisenberg model has an effective low energy description in terms of a spin 3- vector field 



(f) constrained to be of fixed length (j) ■ (j) = constant, with lagrangian d^cj) ■ 
The field is bilinear in the fundamental electron fields, and in our model corresponds to 
= x^(Tx+/-\/2. The constraint on (j) can be imposed by the constraint x~X^ = constant 
since 

0-0= (2) 
Imposing these constraints one finds that 

— * — * 

d^(f) ■ d^(f) oc dfj_x~dnX + irrelevant operators (3) 

which justifies the kinetic term in our model. One can then relax the constraint on x~X~^^ 
and replacing it with a "soft constraint" by including a quartic interaction, as is done for 
the non- linear 0(3) sigma model in two spatial dimensions. 

Another motivation for the kinetic term in our model (at m = 0) was given in 4] based 
on the linear dispersion relation E]f^ = |k| one obtains when expanding around a circular 
Fermi surface. However this involved a modification of, or at best a crude approximation 
to, the density of states. 

Canonical quantization and also a path-integral formulation follow from the euclidean 
action: 

S = d^:sidt {df,Xad,,xt + ^^XaXa) - ^-^^ 9 Xj xt Xj xt (4) 

where df^d^ = + VpV"^. The velocity vp plays the role of the speed of light, and it was 



proposed in [4] that it be identified with the universal nodal Fermi velocity 13j . In the 
sequel we set vp = h = ks = 1 except where indicated. Note that the fields x are treated as 
Lorentz scalars and spin is simply a fiavor. However it is possible to treat the fields as Dirac 



spinors and thereby achieve complete Lorentz invariance|14j]. The quartic term is unique up 
to the sign of the coupling by fermionic statistics, and positive g corresponds to repulsive 
interactions. 

A consequence of the fermionic statistics is that the model has a hidden 5*0(5) symmetry. 
As explained in the Introduction, the appearance of this 5*0(5) is quite different from the 
5*0(5) symmetry proposed by Zhang. This symmetry is manifest if one considers an A^- 
component version with fields Xa^ — which has Sp{2N) symmetry (hence the 

terminology "symplectic fermions"). For N = 2, since there are 4 fermionic fields and 
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consequently a unique 4-fermion interaction, the theory automatically has Sp{A) — SO (5) 
symmetry. 

The S0{5) contains SO{3) and U{1) subgroups which commute and can be identified 
with spin and electric charge respectively. The conserved electric current then corresponds 
to 

Jl^-^Y^ (Xa^MXa + Xa ^MXa ) (5) 

a 

and the fields have electric charge Qf. = ±1. 

The important order parameters for the study of spontaneous symmetry breaking are 
composite bilinears in the fermions. The 4 fields Xa transform under the spinor representa- 
tion of SO{5). The bilinears can be decomposed as4(8)4 = 1©5®10 where 1 is the singlet, 
5 the vector representation, and 10 the adjoint. The singlet is the field J2a XaXt = X~X^ 
and corresponds to the mass term in the action. The 5-vector of fields corresponds to 

1 



^ = 0^ 0e ) = (^X 5^X^, XW' ) 



(6) 



where a are Pauli matrices. The triplet of fields are electrically neutral and transform as 
a spin vector under the SO (3) and serve as magnetic order parameters. The fields (pf on 
the other hand are spin singlets but carry electric charge ±2 and are thus Cooper pair fields 
for superconducting order. The SO (5) invariant product is 

$.$ = 0.0-20+0; (7) 

and the interaction can be expressed in the manifestly SO (5) invariant manner: 

Ant = ^5^-^ (8) 

The momentum expansion of the free fields is 

where a;k = Vk^ + rn? and k ■ x = a;ki — k • x. The canonical quantization of the theory 
based on the lagrangian leads to the canonical anti-commutations in momentum space: 

{ak, al} = {6k, bl} = (27r)2(5(k - k') (10) 



and the hamiltonian is 

— J2 (4,a«k,a + blabk,a) (H) 

The a and b particles have opposite electric charge and can thus be thought of as particles 
and holes. 

A distinguishing feature of our model, in contrast to quartic interactions of Dirac fermions 
for instance, is that the quartic interaction is relevant in the renormalization group (RG) 
sense: the field x has classical mass dimension 1/2 so that the quartic interaction has 
dimension 2 and the coupling g dimension 1. This means that the interactions are important 
at low energies and can lead to non- Fermi liquid behavior. In fact, the model has a low energy 
RG fixed point, i.e. a quantum critical point at g ^ 1/8. This critical point can formally be 
understood as an analytic continuation of the familiar Wilson-Fisher fixed point of the 0{N) 
models, where N = —4, and some critical exponents, such as the anomalous dimension of 
the field can be computed by specializing the known epsilon-expansion results for the 
0{N) model to = — 4{3]. However, many of the important operators such as the order 
parameters $ are composite fermion bilinears, which changes their structure compared to 
the magnetic order parameters of the 0{N) models. Furthermore, $ includes bilinear order 
parameters of charge ±2 that have no counterpart in 0{N) physics. These RG properties 
are summarized in the next section. 

Another interesting feature of our model is that whereas the fundamental interaction 
is repulsive, when one incorporates the momentum-dependent scattering at second order 
in perturbation theory (1-loop) there is an instability toward the formation of a d-wave 
superconducting ground state. This was studied in [j] by deriving a d-wave gap equation 
for momentum-dependent vacuum expectation values for the Cooper-pair fields (pf. This is 
reviewed briefly in section IX where we derive and study the finite temperature version of 
the d-wave gap equation. 



III. UNITARITY OF THE MODEL AS AN EFFECTIVE LOW ENERGY THE- 
ORY. 

The above free hamiltonian (fTT!) in momentum space is obviously hermitian and defines a 
unitary theory, in spite of the fact that the Klein-Gordon type of action normally associated 
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with free bosons was quantized with the "wrong" statistics, i.e. as a fermion. There are no 
negative norm states in the Fock space of a, b particles since there are no unwanted minus 
signs in eq. fllOl) . In Id the free massless symplectic fermion model is normally considered 
a non-unitary c = —2 conformal field theory. The resolution of this apparent contradiction 
can be found in |171], where it was shown that the symplectic fermion can be mapped onto 



the c = 1 unitary theory [l7l| of Dirac fermions at the level of detailed conformal partition 
functions. In light of the above results, the latter fact is not surprising since in momentum 
space the hamiltonians of symplectic and Dirac fermions are identically eq. (fTTj) . 

The issue of unitarity was further addressed in j^, 4], where the main concern is the 
consistency of the interacting theory. Because of the extra minus sign in (Q, the fields 
are not hermitian conjugates of each other, but rather 

X+ = C{x-yC (12) 

where the unitary operator C simply distinguishes particles and holes: CaC = a and CbC = 
—b, with C^C = = 1. The operator C is easily constructed: 

^ = ^"p(^"/(^?^i'>-) 

The interacting theory is thus pseudo-hermitian: = CHC. The free theory is actually 
hermitian in momentum space since it is quadratic in 6's. As before there are no negative 
norm states in the Fock space. 

A pseudo-hermitian hamiltonian has real energy eigenvalues, and can still lead to a uni- 

fi n n n 

tary time evolution, as discussed in [3|, U, Il5|, [IGj. Let us make the following additional 
remarks concerning this issue in this specific context. A unitary time evolution follows if 
one modifies the definition of hermitian conjugation. Define the C-hermitian conjugate as 
A^" = CA^C. Then the S-matrix is unitary with respect to this conjugation: S'^'^S = 1 since 
iJ^'^ = H. However the C-conjugate hermitian conjugation implies a modified definition of 
the inner product: (ip'lip) — > {tp'\ip)c = (V'' I C*| ■?/'). Under this new inner product, states with 
an odd number of b particles now have negative norm. A completely equivalent description 
that dispenses altogether with the operator C amounts to removing the minus sign in the 
expansion of gQ- &■ This gives a hermitian hamiltonian since now x~ = (x^)^- 

However this also gives the modified relation {6k, b]^,} = — (27r)^(5(k — k') and this leads to 
the same conclusion that the b particle sector contains negative norm states. 
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We now argue that this model can still yield a consistent theory at low energy. Let us 
use the description based on the C-operator; the same arguments apply to the manifestly 
hermitian description without it. The probability that an initial state \i) evolves to a final 
state I/) is given by 

Pu - jfl^ (14) 

(^|^)c(/|/)c 

States can be classified according to their C eigenvalue of ±1, where C — —1 corresponds 
to negative norm. One sees from the above formula that the probability Pif is positive so 
long as the states \i),\f) are both of either positive or negative norm. The problematic 
transitions with negative probability thus arise if the matrix element {i\S\f) is non-zero for 
states of mixed norm. However at low energies compared to the mass m, these transitions are 
forbidden as a result of energy and charge conservation. For example, suppose the initial 
state consists of only a-particles. Because of charge conservation, processes that change 
the number of b particles necessarily involve creation of a, b pairs; for instance aa — > ab is 
not allowed. The process a aab is consistent with charge conservation, but not allowed 
kinematically, so one needs to consider at least aa aaab. Energy conservation shows 
that the energy of the initial a particles must be at least 2m in order for the process to be 
kinematically allowed. Thus at low energies compared to the mass m no 6-particles will be 
produced from a state containing only a particles. A variation of this viewpoint interchanges 
the role of b and b\ so that particles are created by the 6's and have negative energy. One can 
now imagine that all the negative energy states are filled according to the Pauli principle. 
Since there is a gap of 2m between the negative and positive energy states, at low energies 
compared to m no ab pairs will be created out of the Dirac sea and one can deal only with 
the a-particles. 

As we describe in the sequel, the mass m is on the order of the cut-off and can be quite 
large, typically greater than lOOOX for HTSC materials, so the low energy approximation 
is expected to be good for a large range of temperatures. At very high temperatures, our 
model necessarily breaks down and must be replaced by another cfTcctivc theory that restores 
unitarity. For instance, the cross-over to ordinary Fermi liquid behavior in for example the 
specific heat is not built into the model. 

Another way to obtain a low energy unitary theory involves the introduction of a chemical 
potential, which is necessary at finite density. Particle number is not conserved here; instead 
a chemical potential couples to the electric charge Qe, i-e. H ^ H + jiQ^. The free 
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hamiltonian is then 



H ^ 




(15) 



Because of the difference in electric charge between the a and b particles, at a given k there 
is a gap of 2/1 between the energies of the a's and 6's. Let us suppose that /i is chosen so 
that all of the 6-states are filled. Typically will be comparable to the Fermi energy and 
typically very large. If /i is large enough, then because of the energy gap 2/1 at low energies 
we can effectively deal only with the a-particle sector. Alternatively, one can view the a 
sector as an in-acccssiblc high energy sector when fi is large enough, and deal only with the 
low-energy b sector. Either way, as discussed above, all probabilities Pif will be positive. 

IV. RENORMALIZATION GROUP PRESCRIPTIONS 

In all approaches to the renormalization group there are two cut-offs, a fixed upper cut- 
off Ac typically related to the inverse lattice spacing, and a lower cut-off A < Ac which 
is the running RG scale. Renormalization effectively removes the degrees of freedom with 
energy between A and A^. In the high-energy physics context, the cut-off Ac is unknown, 
so one normally performs the renormalization by sending Ac to infinity and incorporating 
counter-terms to cancel ultra-violet divergences. In D = 4 spacetime dimensions this can be 
conveniently and systematically performed by using dimensional regularization in D = 4 — e 
dimensions, which leads to the epsilon-expansion technique for studying models in D — 3. 
In the present context, since the cut-off Ac is in principle known and physical quantities 
depend on it, it is more appropriate to perform the RG directly in D = 3 with explicit 
cut-offs rather than use the epsilon expansion. 

The dependence of the couplings on the running RG scale A can be determined by 
considering an upper cut-off A only. Consider first the coupling constant g. The 1-loop 
correction to the vertex is cancelled by a shift of g, i.e. by defining g{A) — g + 5g where 




(16) 



When A > m, g(A) ^ g- 8g'^/A. Thus 



Adj,g{A) 



8g' 



(17) 



A 
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Let us define 

g = Kg (18) 
Then the beta function for the dimensionless coupling is 

A5a?=-?+8?2 (19) 

The above beta-function eq. f|T9|) has a low energy fixed point at ^ = 1/8. It is known 
that this fixed point survives, with small corrections, to higher orders jsl- We can now 
integrate the RG flow and incorporate the initial data at the high energy cut-off A^. It will 
be convenient to introduce the dimensionless variable x: 

X = 1/g (20) 

The fixed point occurs at x^: = 8. Let be the value of the coupling at the scale A^, and 
Xo = 1/^0- Then the solution to the RG flow equation f[T9l) takes the simple linear form in 
x: 

^ = (x. = 8) (21) 

iV(i 

There are two cases to consider depending on whether the coupling is strong [xq < x*) or 
weak {xq > x*) at short distances. Based on the relation between x and hole doping described 
in section VI, we will refer to Xq < x < x,,, as the underdoped region and x* < x < Xq as 
overdoped; in both cases A > 0. Furthermore, in order to clearly distinguish the two cases 
we will refer to xq in the overdoped region as xq- (What we refer to as the overdoped region 
is close but not identical to the usual terminology; the latter refers to the region beyond the 
maximum Tc.) 

It is clear that our model is characterized by a single fundamental energy scale set by the 
cutoff Ac. Since A^ has units of a wave-vector k, this energy scale is 

Eo = hvpAc = ksTo (22) 

where Ac should be proportional to the inverse lattice spacing. 

Before turning to mass renormalization, we now determine the 1-loop corrections to the 
anomalous dimension of the order parameters which are composite fields. Consider first the 
singlet operator x^X^- To first order in perturbation theory, 

dH 1 



{x-X^{0)---) = 8n'g 



(27r)3 (£2 + m 



2^2 



(X-X^(O)---) (23) 
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This contribution is cancelled by renormalizing the operator as follows: 

X'X^ ^ Z^-^+X'X^ (24) 

where ^x"x+ = 1 + ^^x"x+ ^^"^ 

SZ,-,. = -8n^gj (25) 
The anomalous dimension 7x-x+ is then 

7x-x+=A9AlogVx+ = -4? (26) 

where again we have assumed A ^ m. Repeating the above calculation for the S0{5) vector 
of order parameters one finds 

75 = 4? (27) 

Let [O] denote the scaling dimension of O at the fixed point, including the classical 
contribution of 1/2 for each x field. Using = 1/8 at the fixed point, one obtains to 1-loop 

Ix-Xn ^ \, m = IX^^X^} = IxUt^ = IxlXjl ^ I (28) 

Thus one sees that the singlet is shifted down from the classical value, whereas the 5- vector of 
order parameters is shifted up. The above values agree exactly with the 1-loop approximation 
in the epsilon-expansion obtained in ^ . We point out that although the actual values of the 
fixed point value and the functions 7(5') differ in the calculation performed here compared 
to the epsilon expansion, the results for the anomalous dimensions at the critical point agree 



to 1-loop. From the two-loop results obtained in [s], one sees that higher order corrections 
are reasonably small: changes to 5/8 and [$] actually remains at 3/2 at the next 

order. 

Now we consider mass renormalization, which will be important in the sequel. The 1- 
loop correction to the self-energy shown in Figure [2] is a negative correction to m^. This is 
cancelled by —>■ m^(A) = + dm? where 

5m^ = 871^ g ^ /'f, ^ , (29) 
(27r)3£2 + ^2 ^ ) 

The term on the RHS is extracted by taking a derivative with respect to and then 
setting = 0: 

5m^ = Agm^/A (30) 
14 



Defining 

m = Am (31) 

one finds 

Ad^fh^ = (-2 - 4^)m^ (32) 

The — 2m^ term simply corresponds to classical dimension 2 of and 4'g the quantum 
correction. Since is the coupling for x~X~^ this implies that the anomalous dimension of 
X^X^ is —4:'g, in agreement with eq. (!26|) . 




FIG. 2: One- loop correction to the self-energy. 



V. DYNAMICAL MASS GENERATION AND ITS IDENTIFICATION WITH 
THE PSEUDOGAP 

Although the mass m is classically zero in the connection to the 0(3) sigma model 
description of the Heisenberg AF described in section II, the operator x~X~^ is relevant and 
not forbidden by any symmetries so that in general it will be generated when one includes 
interactions. In this section we determine this dynamically generated mass to lowest order 
and propose that it be identified with the HTSC pseudogap. 

If the original theory is massless, at 1-loop the correction to coming from the diagram 
in Figure [2] is negative and equal to —871"^ g J -^ppji- Since the propagator becomes 
l/(^^ + m^) when m 7^ 0, a self-consistent equation for m is the following 

m = -Stt g / . rr (33) 

J (27r)3 + ^ ^ 



The one point function for the singlet is then 



d^£ 1 m 



2 



where we have used the equation (1331) . The two above equations can also be derived from a 
Hubbard-Stratanovich trans ' 
to (x X^)^) ^-s was done in 



brmation based on the fact that the interaction is proportional 
12l- 



15 



In order to make sense of eq. (133!) and obtain solutions, the mass renormalization dis- 
cussed in section IV needs to be taken into account. Taking the limits of integration to be 
from zero to A, eq. becomes 



m 



-AgK + Agmtaxr'^ Kim (35) 



From eq. (129|) with m = on the right hand side, one obtains 5m^ = AgA, therefore the first 
term on the RHS above can be absorbed into m^, and this is consistent with renormalization. 
The equation for m now becomes 

m = 4^tan~^l/m (36) 

and has real solutions which are easily found numerically. 

There are two useful analytic limits to the solution of eq. (!36|) . When ^ is large, m is 
also large and the solution is approximately: 

A 



m 



2V^, (? large) (37) 

When 'g is small, rh is also small and tan~^ 1/m ^ 7r/2. Thus in this limit one has 

^ 2Tcg, {x'X^) ^9 (? small) (38) 



It is interesting to note that if one sends the cut-off to infinity in eq. fl33|) and performs the 
integral by analytic continuation in the spacetime dimension, then one obtains the solution 
( !38l) . This means that the analog of dimensional regularization here overestimates m, since 
by eq. (1361) . m < 271 g. In the underdoped region the approximation (1371) is considerably 
better than ( l38l) and we will use it in places below. 

The mass m corresponds to an x-dependent energy scale 

Epg{x) = hvpTn^x) = kBT*{x) (39) 

where as before x = 1/^ and m{x) is the solution to eq. fl36|l . In units of the fundamental 
scale £^0- 

= — = m{x) — = m{x)— — — (40) 
£"0 To Ac X* - xq 

In the next section we will relate x to doping, thus the energy scale T* is doping dependent. 
It should be emphasized that T* simply corresponds to the temperature independent energy 
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scale Epg and is thus not a real temperature; however as we will show in the sequel it can 
correspond to a cross-over scale in the real temperature. 

A non-zero mass m clearly corresponds to a gap in the density of states since the 1- 
particle energies are -E^ = v^kM-^r?. We discuss this further at the end of this section. 
As a model of HTSC we thus identify the mass m with the pseudogap energy scale, i.e. 
m = Epg = T*. As we will show below, the thermodynamic properties of our model also 
support the identification of T* with the HTSC pseudogap. 

A plot of T* verses x is shown in Figure [31 Note that T* is close to linear near the critical 

point at X* = 8. Furthermore, the pseudogap is smaller on the overdoped side. We point 

out that the re-appearance of the pseudogap on the overly doped side is contrary to what is 

normally observed. 

T*/To 
1.0 r 




FIG. 3: The pseudogap T* as a function of x for xq = 4 and xq = 16. 

Since the overall energy scale Eq depends both on vp and A^ in attempting to compare 
with the HTSC data it is more useful to fix Eq by using the pseudogap. It is known 
experimentally that at zero doping in the region of the AF phase Epg is approximately the 
AF exchange energy J. Since zero doping occurs at xq = (see the next section), we can 
identify Epg{xo) = m(xo)-E'o ~ J- Since m(xo) is of order 1, Eq ~ J. For the cuprates 
J/kB ~ 1300 - UOOK. 

An estimate of the density of states p{E) as a function of energy E can be obtained as 
follows. One has 

Pik)dk = -^,=^ = PiE)dE (41) 
17 



If we approximate k in the above equation as the Fermi wave- vector kp, then 

1 



kp [ dE\ ^ kp 



271 \dk J 271 ^1 -m2/E2 

This behavior is quahtatively shown in Figure HI 

p{E) 



kp 
2-K 




E 



(42) 



FIG. 4: Density of states as a function of energy E. The dotted line represents a sketch of what is 
measured using ARPES. 

Angle-resolved photoemission (ARPES) measurements show a clear depression in the 
density of states and the leading edge of this depression is characterized by the pseudogap 
scale T*. This behavior is shown schematically in Figure H] as a dotted line. Thus one 
qualitative difference with what is observed is that the gap in our model is a clean one. 
Another difference is that due to the rotational invariance of our model, the pseudogap 
is isotropic, whereas the observed pseudogap shows lattice effects. It is anisotropic in the 
direction relative to the Fermi surface: the pseudogap is largest in the anti-nodal (0, tt) 
direction and vanishes along "Fermi arcs" in the nodal (vr, vr) direction. 



VI. HOLE DOPING 



Since renormalization removes degrees of freedom between A and Ac, this suggests that 
doping can be varied by varying the cutoff A. In ^ it was argued that a measure of 
hole doping can be defined based on the 1-point function Since this one point 

function is proportional to A by dimensional analysis, we are led to identify doping p with 
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the dimensionless quantity 



P=^((x"X^)a.-(x"X+)a) (43) 



for some constant c. In the underdoped regime, the 1-point function is approximately given 
byeq. dSH): 

p(-)-4fi-r)-4(^) (44) 



where we have used eq. (1211) . Half-fiUing then corresponds to A = Ac, i.e. x = xq. Thus, in 
the approximation we have made, plots of various physical properties as a function of doping 
p is simply related to plots as a function of the inverse coupling x by rescaling and shift of 
the X-axis. We choose to plot against x since this more clearly reveals the RG properties. 

We can give a rough estimate of the constant c following an argument made in {4]. 
Consider a lattice fermion model where S'x is the local spin variable at lattice site x, which 
is bilinear in the fermion operators. One has 

^ ^ 3 3 

'5-'5 = --n|n| + -(nt + n|) (45) 

where ?t,| j is the number of fermions of spin up or down at each site. From this relation 
one sees that at half-filling n^ni = 0, rij + n| = 1 and 5*^ = 3/4, i.e. 5* is constrained to 
be a spin ^ vector and the model can be mapped onto the Heisenberg model. Dividing the 
above equation by the volume squared and taking the infinite volume limit, one finds that 
the right hand side is — fpfP; where p^^i are number densities. 

In our continuum model, S is represented by the bilinear and one has the identity: 

^■^=-\{x-X^f (46) 

Comparing with the continuum limit of eq. fj45l) . one identifies = (Pt +Pi)/v^5 which 

corresponds to c = \pi. At the critical point this gives pcrit ~ -14. Since this is only a rough 
estimate, one can alternatively fix c in principle by fitting to experimental data. 



VII. SPECIFIC HEAT 

An approximation to the specific heat can be made based on including just the effects of 
the dynamically generated mass. A gas of 4 types of fermionic particles with single particle 
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energies c^k = v^k^ + m"^ has the free energy per volume at temperature T: 

•^=-4Ty[|'^^log(l + e--/^) (47) 

One easily sees that T j IS?^ is a function of the two dimensionless variables T/Tq and T/T* 
where as in section V we have identified T* = m. 

At very low temperatures T ^ To we can approximate the log as e~'-^'^l'^ and effectively 
send the cut-off Ac to infinity, and T jT'^ becomes a scaling function of T* /T . The result is 

..-?i!.-(..i:) 



The entropy density is then: 




(49) 



and the specific heat: 



C=-r^«Hl!e--'T 1 + 11 + Ifn +i fill 1 (50) 



TT 



It is convenient to define the quantity = C jT since for a Fermi liquid 7 is a constant. The 
primary feature of our model is that 7 oc T at very low temperatures compared to T* and 
To, which is ultimately attributed to the relativistic nature of the model. 

At higher temperatures there is a crossover to a different behavior. Since Ac can be scaled 
out of the eq. (H7|) . the cross-over temperature is the pseudogap temperature T*. Below we 
plot the entropy and 7 as a function T for various "doping" x. The plots are in terms of the 
dimensionless quantities T/Tq, S'/A^ and 7/ Ac. One clearly sees the crossover at T = T*. 

Experimental data for 7 is shown in Figure [71 In comparing with our results, it is 
important to bear in mind that the data contains contributions from the quasi-particle 
excitations in the SC phase, which explains the peaks to the left, and such effects are not 
included in our calculation. Thus, our curves should be compared with the data to the 
right of the SC peaks. At temperatures below To one sees a reasonably good qualitative 
agreement with the behavior computed above: for T <^ T*, 7 oc T, with a crossover to a 
different behavior at T*. The dependence on doping is also qualitatively correct, i.e. the 
peaks of the curves move toward the left with increased doping. One difference is that 
whereas our computed 7 goes to zero at high temperature, the experiments indicate that 
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FIG. 6: 7 = C/T as a function of temperature at various "doping" x for xq = 4. The larger 
amplitude curves correspond to larger doping. The vertical axis is the dimensionless quantity 
7/^0- 

at high temperatures 7 approaches the Fermi-hquid resuh, i.e. 7 approaches a constant. 
Our model of course does not crossover to a Fermi hquid at high temperatures, because as 
explained in section III it is expected to break down at high enough temperatures, so at 
best it may describe temperatures up to and above T* . 
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T(K) 



FIG. 7: Experimental data for 7(T) for Yo.8Cao.2Ba2Cu306+x (from 
VIII. NON-ZERO MAGNETIC FIELD 

In this section we compute the spin response to a magnetic field and also the magnetic 
field dependence of the specific heat in the same approximation we made in the last section, 
i.e. we only consider the effects of the dynamically generated mass. Both these quantities can 
be studied by adding a term to the action f (Pxdt ^J2h ■ cj). Comparing with the calculation 
of the effective potential = in [4|| , one sees that the free energy at zero temperature is 

where as before, cj^ = + m?. 
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A. Spin response 



The one point function (0) in the presence of h can be obtained from the first derivative 

of the logarithm of the partition function Z = exp{—J-'V/T) where V is the volume. For a 

single component of 0, one obtains 

_]_d£ _ J_ f dujd^^ / 1 1 \ 

^"^^ ~ V2dh ~ ,/2j {2Try [uj^ + ujl-h uj^ + ool + h) 

At finite temperature T, uj becomes a quantized Matsubara frequency uj = 2'kvT where 

V & Z + 1/2 and one makes the replacement: 



2n 



We need the identity 



This gives 



T y — ^ = ltanhf— ) (54) 

^ a^ + (27iuTY 2a \2T ) ^ ' 

id)) = — - / - — — tanh ^ (h -h) (55) 

2^J {2-KY\^ui^ 2T ^ 'J ^ ^ 

The linear response for small h follows from Taylor expanding the integrand: 

{$) = Xs{m,T)h (56) 

where 

y=(m, T) = — ^— / - — --^(tanh^^H tanh — ^ — 1 ) (57) 

Plots of function of temperature are shown in Figure [HI One sees that the spin 

response is quenched when T < T*, consistent with experiments. 

B. Specific lieat 

Define 

Ajr = jr(/i) _ jr(o) (58) 
Integrating the identity eq. ( l54l) one can show 



J 

(59) 
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FIG. 8: The spin response as a function of temperature at various doping variable x for xq = 4. 
The vertical axis corresponds to the dimensionless quantity Xs ' Eq. 

Therefore up to a constant that is independent of T and does not affect the specific heat: 



= y [ — (i+e-^/^)^ — ) 

One can then define 

^7 = (61) 
Plots of A7 as a function of T for different x, /i are shown in Figures O [101 (These plots 
are in arbitrary units since we have not specified the strength of the magnetic field.) At 
very small h, A7 oc /i^. 



IX. FINITE TEMPERATURE D-WAVE GAP EQUATION AND THE EFFECT 
OF THE PSEUDOGAP. 

A. Zero temperature 

It was shown in 4| that the 1-loop corrections to scattering of Cooper pairs leads to a 
d-wave superconducting gap of the form 

g(k) = 6\kl - kl) = 6^k^ cos(20) (62) 

where g(k) is a Fourier transform of Cooper pairing order parameters (x^(xi)x^(x2)). Be- 
fore turning to the main topic, some remarks on the d-wave property of the gap are called 
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A7(T) 




FIG. 10: A'y{h,T) as a function of temperature for various h = 0.05,0.1,0.15,0.2 (increasing h 
corresponds to increasing amplitudes on the plot) at fixed (x,xo) = (5,4). 

for. Since the d-wave gap equation derived in {4] was based on a rotationally invariant 
hamiltonian, the most general solution involves a rotation by an arbitrary angle 6q, i.e. 
cos 20 — > cos2(^^ — ^o) in the above equation. On the other hand in HTSC the lattice breaks 
the rotational symmetry and the d-wave gap is oriented with respect to the Fermi surface 
in a specific way, namely the gap vanishes in the nodal direction {kx,ky) = (7r,7r), which 
corresponds to Oq = 0. This can be reconciled as follows. Free particles on the lattice have 
dispersion relation = — 2(cos/c2, + cos ky). Thus the Fermi surface has the symmetry 
kx — >• —kx- Requiring this symmetry fixes Oq = 0. 
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In this sub-section we study the effect of the non-zero mass m on the zero temperature 
d-wave superconducting gap equation. The effect of the mass term is simply the shift 
a;^ — > a;^ -|- and the d-wave gap equation derived in j4| becomes: 

6' = 2g, [ dude ( 1 - u' + k' + rr l \ 

where A;^ = k^. The coupling g2 comes from the 1-loop scattering of Cooper pairs and is 
given by 

^ SttV f dH 1 

For simplicity we incorporate the cut-off as follows: \uj\ < oo and k'^ < which is more 
appropriate for comparison with the finite temperature version we consider below. 

Re-expressing the gap equation in terms of dimensionless quantities by rescaling k — >■ Acfc 
and oj Accu, one obtains 

1 - , + + \ (65) 

where m' = mj K^. and ^2 = 5'2A|?. 

Since eq. ( IMl) is ultra-violet convergent, it can be approximated by letting the upper 
cut-off go to infinity, giving g2 = ng'^/AOm^. Incorporating the RG prescriptions of section 
IV, g = Kg and m = Am, one finds that the parameters in the gap equation (1^ are the 
following: 

92 = — , m = —m (dd) 
^ 40 m5 V A / Ac 

where the ratio A/ Ac is given in eq. ( 1211) . Finally, since under a RG transformation 6k — >■ 
6Ak/Ac = the physical gap in the theory at RG scale A is 5' = A5/ A^. The x-dependence 
of the solutions 5 arises from the x-dependence of ^2 and m' . A plot of 5' as a function of 
X is shown in Figure [TTl It is important to point out that the interpretation of the critical 
point at X* presented here differs from the original proposal in ^|, in that in the present 
work we extend x beyond x,, by introducing xq, and this places the maximum value of the 
gap near the critical point x* = 8. (We have effectively patched together what was referred 
to as Type A and B in [4].) Further justification for this location of the critical point is 
based on the calculation of T^, which reaches a maximum value near the critical point, as 
explained below. 
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FIG. 11: Plot of (x) vs x for xq = 4 and xq = 12. 

One can easily verify numerically that the pseudogap competes with SC in the following 
sense: if one artificially increases the mass m, then the value of the gap 5 decreases, and 
eventually vanishes for m too large. (This was shown already in j^.) 



B. Finite temperature 

In order to derive the finite temperature version of the above d-wave gap equation, we 
start with the un-integrated form derived in ^: 

[ dujdPk! g(k') 

= w'^"'- * i.^^-+^if + ,[\^y- '''' 

where the kernel is 

^(k, k') = -8TT^g2k^k'^ cos 2(0 - 9') (68) 

The equation fl63l) is obtained upon performing the angular integral. 

At finite temperature, u becomes a quantized Matsubara frequency = 2t:Tu, where T 
is the temperature and z/ is a half-integer, i.e. v E Z + 1/2. As before, the integral j duj/2'K 
is replaced with T^^. One needs the identity: 

tY— \- = ^ImV -— = -Im ( ^tanh f^^)^ (69) 

where uj^^q = — iq. 



27 



Due to the specific form of the kernel G, the solution to the equation (1671) is of the d-wave 
form fl62l) (up to an arbitrary rotation) where 6 satisfies the integral equation 



= g2 I dkdd k^cosi2d)lm ( ^tanh ( 

J V^k,5 V 



where 



^ i^'^k'^ cos 26 



(70) 



(71) 



Finally, the physical temperature at the scale A follows from the RG transformation T — »• 
TA/A,. 

Solutions of the above equation are 5(x,T). One can easily verify numerically that as 
the temperature goes to zero, one recovers the solution 5{x) to the zero temperature gap 
equation fl65|) . One also finds that as the temperature is raised there are no solutions to the 
above equation for T > and this defines the x-dependent critical temperature Tc. This is 
shown in Figure lT2l 

5'{T) 




0.01 0.02 0.03 0.04 0.05 0.06 0.07 



T/To 



FIG. 12: The d-wave gap 5'{T) as a function of T for x = 7.5, xq = 4. 

The critical temperature Tc as a function of x on both sides of are shown in Figure 
lT3l It turns out that Tc at the critical value x = x^, i.e. T* = Tc(x=k), is universal in that it 
only depends on the overall scale Tq and not on Xo,Xq since the scale factor A/Ac vanishes 
at this point. Using the estimated relation between x and doping in section VI, the critical 
point occurs at doping Pcrit ~ -14. The dome shape of Tc is a property of the mathematical 
structure of the gap equation and nothing universal is happening at the termination points 
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of SC on either side of the critical point. Numerically we find that at the critical point 



-f ^ .084, 
To 



To^'fx* 



.268 



(72) 



We also find numerically that the maximum value of Tc occurs close to the critical point so 
that T^^^ ^ .084To, i.e. T™^^ is simply proportional to fundamental energy scale Eq. 

As argued in section V, To should be identified with the anti-ferromagnetic exchange en- 
ergy J at half-filling. For To = ISSOi^', this gives T™**^ WZK, which is quite reasonable for 
HTSC It should be emphasized that the above T^ is intrinsic to the two spatial dimensions, 
i.e. does not involve any kind of inter-planar energy scales. 

T*/To, Tc/To 




FIG. 13: Plot of Tc and T^, vs x for xo = 4 and xq = 12. 



X. CONCLUDING REMARKS 

In summary, we have further developed the interacting symplectic fermion model in two 
spatial dimensions by studying a dynamically generated relativistic mass and by including 
a finite temperature. This allowed us to study some fundamental properties of the model, 
such as the specific heat and spin response, which clearly show non-Fermi liquid properties. 
As a simplified model of HTSC, we identified the pseudogap energy scale with the zero 
temperature relativistic mass m, and pointed out some close parallels with the observed 
phenomenology of the pseudogap. 

We studied the effects of the pseudogap and finite temperature on the d-wave gap equa- 
tion. In this model, the pseudogap clearly competes with superconductivity as a distinct 
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phenomenon. Our analysis of Tc suggests that the quantum critical point of our model, 
where the pseudogap vanishes, occurs inside the superconducting dome near optimal dop- 
ing. For an antiferromagnetic exchange energy of J ~ ISSOi^T, solutions of the d-wave gap 
equation give a maximum T^. of about 
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